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CN ■ Introduction 

X> 

| ■ The in-matter three-body problem plays an important role in describing a large variety of 

interesting phenomena in many-body systems. For example, in order to understand the formation 
of bound states in heavy ion collisions, three-body calculations are needed to study the 
modification of the binding energy and wave function of a three-nucleon bound state due to 
j> ■ nuclear matter of finite density and temperature [1,2]. Similarly, studies of the binding energy of 
three quarks are of relevance to the understanding of color superconductivity and phase transitions 
in quark matter [3,4]. 

<N 

Although hole contributions make even non-relativistic in-matter descriptions four-dimensional, 
the numerical complications of a four-dimensional approach makes a three-dimensional reduction 
desirable. In the early work of Schuck, Villars, and Ring [5], equal-time Green functions were 
used to obtain a three-dimensional field theoretic description. To derive their equation for the 
zero-temperature equal time three-body wave function, they approximated the effective 
pair-interaction kernels by terms linear in the physical two-body potentials. Since the exact 
expression for the effective pair-interaction kernel involves an infinite series of higher order terms 
as well [see Eq. (fT2l) l. the linear approximation cannot be considered as satisfactory for the strong 
^ ■ coupling case, e.g., when two-body bound states are possible. Despite this, the model of Ref. [5] 
has remained to the present day the state-of-the art formulation and forms the basis of recent 
calculations [1-4,6-8]. 

The goal of the present paper is to formulate three-dimensional equations for the finite 
temperature in-matter three-body problem, that take into account the full infinite series for the 
effective pair-interaction kernel, so that all possible two-body sub-processes allowed by the 
underlying Hamiltonian are retained. 
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In-matter four-dimensional three-body equations 

The interactions of three identical particles at finite density and temperature are described in 
quantum field theory by the Green function Q defined by 



(2n) 4 § 4 (pi+P2 + P3-pi-p2-P3)Q(p'iP2P3>PiP2P3) = / d 4 y\d 4 y 2 d 4 y? ) d 4 x\d 4 x 2 d 4 x? l 



where ^ and are Heisenberg fields, T is the time ordering operator and 
p = e _ P( ff_ W /Tve~^ H ~v N \ is the statistical operator of the grand canonical ensemble [9]. 
Besides being the central quantity for the description of three-body observables, this Green 
function is also needed to calculate the vacuum properties of the system with the help of the 
dressed single particle propagator; for example, in the four-point interaction model, the single 
particle self-energy diagram is completely defined by particle-particle-hole (pph) Green 
function [10]. 

For non-zero temperature, two types of perturbation theory, so-called "imaginary-time" and 
"real-time", have been derived for Eq. © [11]. Here we shall use the real-time formulation of 
perturbation theory in which the number of degrees is doubled [11,12]; this complication, with 
respect to the zero-temperature case, comes from the sum over the complete set of states (trace) in 
Eq. O (a similar discussion on the basis of the imaginary-time formalism will be presented 
elsewhere). For example, the nonrelativistic free one-body propagator is given by [13] 

d f ip) = l \^L_ + ^LJ (2) 

[p°-io + ie p°-(o-ie\ 

where G) = co p = p 2 /2m — /j, /j is the chemical potential, and n, n are 2x2 matrices whose 
elements are functions of the distribution function /((fl) = (e$ a ± where the upper sign (+) 
is for fermions and the lower sign (— ) is for bosons. The n and n satisfy the following relations 
which define their projection properties: n + n = g, ngn = ngn = 0, ngn = n, and ngn = n, 
where g is a 2 x 2 matrix with elements g\\ = 1, g\ 2 = g2i = 0, and g 22 = ±1. Correspondingly, 
elementary vertices have an extra double-valued index for each particle leg. 

For the identical particle case considered here, the field theoretic expression of Eq. (HI 
automatically guarantees the appropriate symmetry of the three-particle Green function Q . 
Moreover, in the doubled degrees of freedom formalism, the matrix Green function G whose first 
diagonal element is Q , is likewise properly symmetric in the case of bosons, and antisymmetric in 
the case of fermions. One can write [14] 




G = G P Q + G VG 

where Go is the product of single particle dressed propagators d{pi), 

Go{PiP 2 P3,PiP2P3) = d(pi)d(p 2 )d(p 3 ) (2n) 4 d(p' 2 -p 2 ) (2n) 4 b(p' 3 -p 3 ), 
Gq is Go summed over permutations of initial or final particle labels, and 



(3) 



(4) 



v= 1 -(v l +v 2 +v 3 ), 



VKl'2'3, 123) = v(fk'jk)d-\i)5(i',i) 



(5) 
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[(ijk) is a cyclic permutation of (123)] is the two-body potential with spectator particle i. Note 
that in this work we follow the literature and neglect three-body forces (the connected part of V). 
Eq. © is a four-dimensional integral equation whose input two-body interactions originate from a 
model second quantized Hamiltonian, e.g. H = £1 iO\a\a\ + £ 12 34v(1234)a{a2a4a3. For later 
convenience we also define the disconnected four-dimensional t matrices 7] where 

Ti = Vi + ^ViG Q T u T i (l / 2'3,123)= t(fk',jk)d-\i)b(i',i). (6) 



Three-body equal time quasi-potential 

Being a four-dimensional integral equation, Eq. © involves relative times (or relative energies) as 
integration variables, and is therefore not very convenient for practical calculations. For this 
reason we implement a three-dimensional reduction of this equation. To do this we follow the 
current literature and effect this reduction by equating times in initial states, and separately, in 
final states. Thus, in the double-degree of freedom formalism, our central quantity is the two-time 
Green function (G) which is obtained from the four-dimensional Green function G by equating 
times as just described. In momentum space, (G) is given by 

i(G) ( £>P i^, PlP2P 3)./^^^MM^ G( ^ >pW3 ) 

(2n) 2 d(p[° +p ' 2 + p ' 3 -E)8(p 1 +P 2 + P°3-E). (7) 

Similar expressions hold for one and two-body Green functions. In order to define the equal time 
quasi-potential, one needs to use the inverse of the equal-time disconnected Green function (Go); 
however, in the many-body case this inverse may not exist. For example, the equal-time 
three-particle free Green function (Gq) projects onto the sub-space projected by [15] 



9i =n(pi)n(p 2 )n(p 3 )+n(pi)n(p 2 )n(p3); (8) 



as a result, (Gq) cannot be inverted in the full space of momenta. To get around this problem we 
introduce modified Green functions 

Go = (G ) + (ggg - )A, G = (G) + (ggg - *C ) A P (9) 

For non-zero A, Go is not singular and can be inverted. In Eq. ©, ggg is a direct product of g's 
(one for each particle). It is important to note that Go is identical to (Go) in the subspace projected 
by 9t (i.e. 9{gggGo = %Cggg{Go) = (Gq)). The operator A is required to be fully disconnected, 
but can otherwise be chosen according to one's own convenience. For the free case, one can write 
down the inverse of Gq explicitly [15]. The three-dimensional quasi-potential V is then defined to 
satisfy the equation 

G = G D + G VG. (10) 
Similar to Eq. ©, the quasi-potential V is expressible as 

V = ~ (Vi +V 2 + %) , VK1'2'3', 123) = v'M'j) (11) 
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where Vi is a pair interaction with particle i as spectator. It can be shown that 



Vt 

2 =G ° 



(Go yGo) + (Go yGo yGo) - (Goy Go)G 1 (Gq^-Gq) + ... 



(12) 



Even though Eq. (IT2b is an infinite series that is very difficult to sum, what enters the three-body 
Faddeev equations is not Vi but the pair-interaction t matrix, %, defined in terms of the 
quasi-potential % by the Lippmann-Schwinger equation 



7i(l'2'3',123) 



t-d(ij) 



(13) 



The task of constructing % thus appears to be especially formidable. That is why most works on 
this subject keep only the linear term in the input two-body interaction [5]. In our case, this would 
mean keeping only the first term of the series in Eq. (fTZl) . This is a practical but unreliable way out 
of the problem. However, as we show below, there is another way of solving this problem which 
gives the exact and complete summation of Eq. (fT2b followed by an exact solution of the 
Lippmann-Schwinger equation, Eq. ( fT3l . 



Exact three-body equal-time disconnected t matrix 

Defining the three and four-dimensional disconnected Green functions G" and G" as 

G" = Go + ^GoWf, G" = G + ^G V ; G?, (14) 
it follows from Eq. © and Eq. (fT3l that 

G" = Go + -Go^Go, G" = Go + -GoT/Go- (15) 
Taking equal times in the four-dimensional versions of Eq. (fT5l) . it is not hard to see that 

G 7}Go = (G 7}Go>. (16) 

Thus, in contrast to the quasi-potential Vi which is related to the four-dimensional potential Vi in a 
very complicated way, the t matrix % corresponding to the quasi-potential, defined by the exact 
solution of Eq. (fT3l . is connected to the four-dimensional t matrix 7] in a very simple way. 

Using the disconnectedness structure given in Eq. © and Eq. ( fT3l . and repeating the above 
argument for two-particle Green functions, one can show that Eq. (fToT) leads to the result 

Goi' i Go=D 0i t i Doi®{di) (17) 

where Do, is the disconnected two-body propagator of particles j and k defined by analogy with 
Go, and t is the t matrix obtained from the two-body quasi-potential v. The symbol ® denotes the 
convolution integral: a ® b(E) = i/2n J°° 00 a(E — z)b(z) dz. 
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